Measurements of Higher-Order Flow Harmonics in Au+Au Collisions 

at ^/i^ = 200 GeV 
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Flow coefficients ii„ for n = 2, 3, 4, cliaracterizing the anisotropic collective flow in Au+Au 
collisions at ^/sNN = 200 GeV, are measured relative to event planes ^l/n, determined at large 
rapidity. We report v„ as a function of transverse momentum and collision centrality, and study the 
correlations among the event planes of different order n. The Vn are well described by hydrodynamic 
models which employ a Glauber Monte Carlo initial state geometry with fluctuations, providing 
additional constraining power on the interplay between initial conditions and the effects of viscosity 
as the system evolves. This new constraint improves precision of the extracted viscosity to entropy 
density ratio rj/s. 

PACS numbers: 25.75.Dw, 25.75.Ld 



The production of particles in heavy ion collisions at 
the Relativistic Heavy Ion Collider is anisotropic in di- 
rections transverse to the beam. For low momentum par- 
ticles (pt S 3 GeV/c), this anisotropy is understood to 
result from hydrodynamically driven flovi' of the Quark- 
Gluon Plasma (QGP) [i-^]. The strength of the flow is 
measured as Fourier coefficients w„ = ^g*n(0-*Rp)^ — 
2,4, .. where (f) is the azimuthal angle of an emitted par- 
ticle around the z axis defined by the beam, usually at 
midrapidity; ^rp is the azimuth of the reaction plane 
defined by the beam direction and the impact vector be- 
tween the colliding nuclei. The brackets denote averag- 
ing over particles and events. The reaction plane is not 
measurable directly a priori, so the Fourier coefficients 
are determined with respect to the estimated participant 
event planes [if. Recent measurements have primarily 
focused on the even-order anisotropics V2 and W4, evalu- 
ated with respect to an event plane '^2, determined from 
the n — 2 correlation. 

The W2(f4) values obtained this way for a broad range 
of and centrality have been used to extract the specific 
viscosity 77/s (the ratio of viscosity 77 to entropy den- 
sity s) of the hot and dense nuclear matter via hydro- 
dynamic model comparisons 0, These model com- 
parisons, which incorporate the dynamic evolution of an 
early-stage strongly-coupled QGP, together with a late- 
stage hadronic gas, give estimates which span the range 

47r-j~l 2. A conjectured lower bound for the specific 

viscosity is 47r^ = 1 Q The rather large uncertainties as- 
sociated with this range of estimates (100%) are currently 
dominated by the uncertainty on initial state anisotropy 
estimates 0, [3 ■ Specifically, the ends of the range are 
given by two equally successful parameter sets. The lower 
bound value is obtained with a standard Glauber Monte 
Carlo model [l^ . [ll| of the initial state which results in 
smaller initial anisotropy and thus needs less viscosity to 



reproduce the measured final state particle anisotropy. 
The higher value An^ ~ 2, corresponds to a larger ini- 
tial anisotropy in the Color Glass Condensate (CGC) 
inspired Monte-Carlo-Kharzeev-Levin-Nardi (MC-KLN) 
model [13, [l^l of the initial state. 

Recently, significant attention has been given to the 
study of the influence of initial geometry fluctuations of 
the initial state anisotropy, which are typically quantified 
by higher-order generalized "eccentricities" e„ [3, 
with the goal of understanding how such fluctuations in- 
duce anisotropic particle emission, characterized by w„ 
(for odd and even n) 



—— oc 1 + > 2vn cosfn 

n—l 



(1) 



where u„ = (cos(n[0 — 5'n])) ,n = 1, 2, 3, ... and the ^I'n 
are the generalized participant event planes at all orders 
for each event. These recent developments suggest that 
measurements of Vn, especially for n = 3, can yield impor- 
tant additional constraints that provide a more precise 
estimate of ^, as well as resolve the correct eccentricity 
model. 

Here we present results for differential measure- 
ments following Eq. [U for Au-|-Au collisions at 
\/sNN— 200 GeV. We first show how the measured event 
planes correlate across large rapidity gaps, and then show 
resulting u„ moments for midrapidity particles relative 
to those planes. We find that the measured Vn mo- 
ments, in conjunction with hydrodynamical model calcu- 
lations 3 13 J indeed provide new constraining power 
for both the initial state and rj/s. 

The results are derived from ^ 3.0 x 10^ Au-|-Au 
events obtained with the PHENIX detector [l^ during 
the 2007 running period. Collision centrality (related 
to impact parameter) and number of participating nu- 
cleons (A'part) determinations were performed with pre- 
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FIG. 1: (color online) Raw correlation strengths (see text) 
of the event planes for various detector combinations as a 
function of collision centrality. The detectors in which the 
event plane is measured are: (a) RXN North, (b) BBC South, 
(c) MFC North, and (d) MFC South. 



viously described methods [l9\. Event planes were de- 
termined using three separate detector systems: Beam- 
Beam Counters (BBC) [2^, Reaction-Plane Detectors 
(RXN) and Muon Piston Calorimeters (MPC). Each 
detector system has a North (South) component to mea- 
sure at forward (backward) rapidity. The absolute pseu- 
dorapidity (77') coverage for these detectors are 3.1 < 
KbcI < 3-9, 1.0 < h^,xNl < 2.8, 3.1 < I 
The PHENIX drift and pad chambers ^ 
for charged particle tracking and momentum reconstruc- 
tion with azimuthal coverage ip = ir in the central region 
(h'l < 0.35). 

To estimate the event plane xE'n in each detector, we 
generalize to all orders n our earlier procedure for event 
plane determination (see [l9| and especially definitions 
in [23]). For each event plane detector we evaluate 
tan(ri$„) ~ ^ sin(rt0i)/ ^ cos(n(/)i) for the 
subevent estimator $„, where the (jji are the azimuths 
of elements in that detector and the weights Wi reflect 
the energy or multiplicity in that element. Corrections 
for detector imperfections were also employed to ensure a 
uniform response. In general, the hit distributions sam- 
ple virtually all momenta. 

To measure w„, the azimuth of each particle is 
correlated with the VPn, and calculated as Vn{'i'n} — 
(cos(n[(/) - <i>^''s]))/Res('J'n) where is the average of 
the <&„ for North and South subevents and where the 
denominator Res(\E'„) represents a resolution factor de- 
scribed in 1231. This factor corrects d„ for the event- 



by-event dispersion of the Its magnitude can be es- 
timated via the two and three sub-events method [l^ 
in which the correlation between from different sub- 
events is measured. The strength of this correlation is 
generally quantified as (cos(n[$;^ — $^])) for sub-events 
A, B, which measures the cosine dispersion of 'i'n- 

Figure [1] shows the centrality dependence of this cor- 
relation strength (cos(j[$;^ — $„])) for sub-event com- 
binations {A, B) involving different event-plane detectors 
with A77' ^ 5 and A77' ^ 7. The raw correlations are pre- 
sented as measured so that the magnitudes shown should 
be considered specific to the PHENIX detectors involved. 
The systematic uncertainties (not shown) for these cor- 
relations are of similar relative size to those for Vn{'^n} 
discussed below. The uncertainties are correlated across 
centrality and n such that the relative size of these event 
plane correlations can be compared. The magnitudes 
for the odd parity quantities (sin(j[<i>J^ — <&„])), which 
should vanish, are found to be consistent with zero for 
all centrality, j, and $ combinations. Figure [1] pan- 
els (a) and (b) show the two sub-events correlations for 
m = n\ (c) and (d) show the two sub-events correlations 
for m ^ n. The negative correlation indicated in (a) 
is due to the well known antisymmetric pseudorapidity 
dependence of sidewards flow vi , as well as momentum 
conservation 2] . Positive sub-event correlations are indi- 
cated in (a) and (b) for 5*2, 3. 4, with sizable magnitudes 
for 5*2.3 and much smaller values for 5*4. 

The sub-event correlations (cos(j[$^ — $^1)) foi" 7^ 
m are also of interest. Fig. 1(c) confirms the expected 
correlation between ^1 and ^2 (due to sidewards flow), 
as well as that between 5*2 and 5*4 [isj. By contrast, 
Fig. 1(d) shows that there is no significant correlation 
observed between ^2 and ^'3. The order j = 6 is chosen 
to account for the n-multiplet of directions (27r/n) of ^2 
and '^z- The absence of this correlation suggests that the 
fluctuations for ^"3 about "^2 are substantial. This is well 
reproduced by Glauber modeling ij, 2^ and therefore 
support an initial state fluctuation origin of '53 and ^3. 
A small correlation between ^'3 and VPi is indicated in 
Fig. 1(d). While such a correlation seems to be at odds 
with the absence of a 5*2 — 5*3 correlation [cf. Fig. 1(d)], 
we note that — \l/3 correlations need not contribute to 
a residual contribution to '^2 — ^3 correlations through 
5*1. That is, could correlate with and \l/2 in exclu- 
sive event classes. Comparisons using the PHENIX Zero 
Degree Calorimeter, which measures the n ~ \ specta- 
tor neutron event plane [l^] at \rj'\ > 6.5 indicate that 
this correlation has significant yy'-antisymmetry. We de- 
fer further investigation of these correlation subtleties to 
future work. 

Fig. [2] shows results for the midrapidity Vni'^n] foi' 
tracks in the central arms as a function of pr for differ- 
ent centralities. They are from the RXN-defined event 
plane analysis, because this detector has the best resolu- 
tion. The systematic uncertainties for these measure- 
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ments were estimated by detailed comparisons of the 
results obtained with the RXN, BBC, and MFC event 
plane detectors and subevent selections. They are ^3%, 
~ 8% and 20% for W2{^2}, wsi^s}, and W4{5'4}, re- 
spectively, for midcentral collisions and increase by a 
few percent for more central and peripheral collisions. 
Through further comparison of the results obtained with 
the RXN, BBC, and MFC event plane detectors, pseudo- 
rapidity dependent nonflow contributions that may influ- 
ence the magnitude of Wn{^Ti}, such as jet correlations, 
were shown |19j to be much less than all other uncertain- 
ties for V2{'^2} and U4{vE'2}. 

The Vni'^n} values shown in Fig.[2]increase with px for 
most of the measured range, and decrease for more cen- 
tral collisions. f2{^'2} and i'4{\E'4} increase as expected 
from central to semi-peripheral collisions, which reflects 
the increase of e„ in peripheral collisions, walk's} ap- 
pears to be much less centrality dependent, with values 
comparable to t;2{^'2} in the most central events. This 
behavior is consistent with Glauber calculations of the 
average fluctuations of the generalized "triangular" ec- 
centricity £3 Hi, [13 . The Fig. H panels (c), (d), (g), and 
(h) show comparisons of U2{^2} and walk's} to results 
from hydrodynamic calculations. The px and centrality 
trends for both V2{'i^2} and waj^'a} are in good agree- 
ment with the hydrodynamic models shown, especially 
at pt below w 1 GeV/c. 

Figure[3]compares the centrality dependence of W2{^2} 
and wsl^a} with several additional calculations, demon- 
strating both the new constraints the data provide and 
also the robustness of hydrodynamics to the details of dif- 
ferent model assumptions for medium evolution. Alver 
et al. [3l use relativistic viscous hydrodynamics in 2+1 
dimensions. Fluctuations are introduced for two differ- 
ent initial conditions. For Glauber initial conditions, the 
energy density distribution in the transverse plane is pro- 
portional to a superposition of struck nucleon and bi- 
nary collision densities; in MC-KLN initial conditions 
the energy density profile is further controlled by the 
dependence of the gluon saturation momentum on the 
transverse position 1^ |l3| . These two models of the ini- 
tial state are paired with two different values of 47r^ = 
1 and 2, respectively. Both values reproduce the mea- 
sured W2{^2} equally well and the viscosity differences 
reflect the different initial £2. The two models have sim- 
ilar £3, and thus the larger viscosity needed in the MC- 
KLN model corresponds to lower V3 than for Glauber. 
Consequently, our measurement of f3{^'3} helps to dis- 
entangle viscosity and initial conditions. The efficacy 
of these 2-1-1 hydrodynamic results for Glauber initial 
conditions are confirmed further calculations with dif- 
ferent model assumptions. Petersen et al. [2^ deter- 
mine a Glauber initial state event- by- event, translat- 
ing through pre-equilibrium with the UrQMD transport 
model [271 |28|, then evolving the medium with ideal 
QGP hydrodynamics (77/s = 0), and finally switching to 
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FIG. 2: (color online) Vn{^-n\ vs. pT measured via the reac- 
tion plane method. The curves a re p redictions from two hy- 
drodynamic models: Alver et al. and Schenke et al. [ITI ]. 



a hadronic cascade (which has an effective viscosity) as 
regions become dilute. B. Schenke et al. 17] use event- 
by-event Glauber initial conditions, evolved with ideal 
3-1-1 dimensional hydrodynamics, which includes the ef- 
fects of viscosity in the plasma phase. 

All of these models are compared with W2{^'2}, and 
■^si^s} data as a function of A^part in two px bins. All 
calculations describe V2{'^2} well at px — 0.75 GeV/c. 
Deviations from hydrodynamics should be expected in 
peripheral collisions, where nonequilibrium effects may 
be large. At higher px, differences between the calcula- 
tions become more apparent. All models still agree with 
V2{'^2}, including MC-KLN initial conditions. How- 
ever, the lower panels of Fig. [3] show the constraining 
power of W3{4'3} and that the calculated results from 
viscous hydrodynamics, with MC-KLN initial conditions 
and 47r^ = 2, lie significantly below the data. This is 
more apparent in the higher px bin, even in the most cen- 
tral collisions. Therefore, our comparisons suggest that 
the combination of MC-KLN initial conditions in concert 
with 47r^ = 2 is disfavored by our new i'3{^'3} measure- 
ments. By contrast, the results from the hydrodynamical 
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FIG. 3: (color online) Comparison of Vn^^n} vs. A'part mea- 
surements and theoretical predictions (see text): "MC-KLN 
+ 47rf = 2" and "Glauber + A-k^ = 1 (1)" T^; "Glauber + 
47rf = 1 (2)" "UrQMD" [Hf;. The dashed lines (black) 

around the data points indicate the size of the systematic 
uncertainty. 

calculations which employ Glauber initial condition fluc- 
tuations and An^ — 1 show relatively good agreement 
with the W2,3{4'2,3} data. The exact statistical signifi- 
cance of these constraints should be determined through 
a global fit procedure, after quantitative accounting of 
the breakdown of hydrodynamics in peripheral collisions, 
as well as the systematics associated with the averaging 
of eccentricity fiuctuations within these models are fully 
understood. From our data it is already clear that the 
higher order moments w„{^ri} for n > 3 provide an im- 
portant avenue for constraining different physical prop- 
erties of the QGP. 

In summary, we have presented participant event plane 
correlations and differential measurements of 
for n = 2, 3, 4 for charged hadrons using the generalized 
event plane method. The higher order harmonic mo- 
ments W3{5'3} and U4{5'4}, as well as strong correlations 
between the higher order event planes across a large ra- 
pidity gap of Ari'>7, 

provide evidence that the initial state has transverse 
geometry fiuctuations of the generalized eccentricities 
which are then propagated in the hydrodynamic evolu- 
tion of the plasma afterwards. This evidence, includes 
(1) a lack of correlation between the measured event 
planes of order n — 2 and 3 as predicted by Glauber 
modeling, assuming correlations of the event planes with 
the generalized eccentricity, (2) proper description of the 
shapes of the pr dependence in the low pr region by 
hydrodynamic calculations, and (3) agreement with sev- 
eral different initial state -I- hydrodynamic models across 
centralities for order w„{5'„}. The combined results for 
V2,3{^'2,3}, in concert with initial hydrodynamic model 
calculations now suggest that the large (100%) uncer- 



tainty previously associated with the extraction of An- 
can be significantly reduced. Within the models consid- 
ered, 477^ ~ 1 is favored, which is close to the conjectured 
lower bound for the specific viscosity. 
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